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Abstract 



Polarized parton distribution functions are determined by using world data 
from the longitudinally polarized deep inelastic scattering experiments. A 
new parametrization of the parton distribution functions is adopted by tak- 
ing into account the positivity and the counting rule. From the fit to the 
asymmetry data A\, the polarized distribution functions of u and d valence 
quarks, sea quarks, and gluon are obtained. The results indicate that the 
quark spin content is AS =0.20 and 0.05 in the leading order (LO) and 
the next-to-leading-order (NLO) MS scheme, respectively. However, if x 
dependence of the sea-quark distribution is fixed at small x by "perturba- 
tive QCD" and Regge theory, it becomes AS = 0.24 ~ 0.28 in the NLO. 
The small- x behavior cannot be uniquely determined by the existing data, 
which indicates the importance of future experiments. From our analysis, 
we propose one set of LO distributions and two sets of NLO ones as the 
longitudinally-polarized parton distribution functions. 



PACS number(s): 13.60.Hb, 13.88.+e 
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I. INTRODUCTION 



For a long time, deep inelastic scattering (DIS) of leptons from the nucleon has served 
as an important tool for studying the nucleon substructure and testing the quantum 
chromodynamics (QCD). Structure functions of the nucleon have been measured with this 
reaction in great precision, which often provides a firm basis of a search for new physics 
in hadron collisions. In addition, basic parameters of QCD such as a s or Aqcd have 
been obtained from the Q 2 dependence of the structure functions. Consequently, hadron- 
related reactions at high energies are described by the parton model and perturbative 
QCD with reasonable precision. 

The measurement of the polarized structure function g\{x, Q 2 ) by the European Muon 
Collaboration (EMC) in 1988 [1] has, however, revealed more profound structure of the 
proton, which is often referred to as 'the proton spin crisis 1 . Their results are interpreted 
as very small quark contribution to the nucleon spin. Then, the rest has to be carried by 
the gluon spin and/or by the angular momenta of quarks and gluons. Another consequence 
from their measurement was that the strange quark is negatively polarized, which was 
not anticipated in a naive quark model. 

The progress in the data precision is remarkable in post-EMC experiments. The final 
results of the Spin Muon Collaboration (SMC) experiment [2] have been reported, and its 
value of A\ at the lowest x has decreased in comparison with their previous one [3]. The 
final results of high-precision A\ and Af data have been presented by the Stanford Linear 
Accelerator Center (SLAC) E143 collaboration [4], and they consist of more than 200 
data points. Moreover, the measurement of g\ (x, Q 2 ) with the pure hydrogen target has 
been carried out by the HERMES collaboration [5]. In addition to such improvements in 
the data precision, new programs are underway or in preparation at SLAC, Brookhaven 
National Laboratory - Relativistic Heavy Ion Collider (BNL-RHIC) [6], European Or- 
ganization for Nuclear Research (CERN) [7], etc., and results are expected to come out 
in the near future. On the other hand, theoretical advances such as the development of 
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the next-to-leading-order (NLO) QCD calculations of polarized splitting functions [8,9] 
stimulated many works on the QCD analysis of polarized parton distribution functions 
(PDFs) [10-16]. There is an attempt to obtain next-to-next-leading order (NNLO) split- 
ting functions [17] and we can expect further progress in the precise analysis of polarized 



In this paper, we present an analysis of world data on the cross section asymmetry A 1 
in the polarized DIS processes for the proton, neutron, and deuteron targets. We formed 
a group called Asymmetry Analysis Collaboration (AAC), and our goal is to determine 
polarized PDFs, Afi(x, Q 2 ), where i = u,d, s, u, d, s, and g. Another possible approach 
is to parametrize structure functions, g^(x, Q 2 ) (N = p, n, and d), which can be expressed 
as linear combinations of the PDFs. In the analysis and predictions of the cross section 
asymmetry in polarized hadron-hadron collisions, however, what we need are polarized 
PDFs rather than structure functions, because the contribution of each quark flavor is 
differently weighted in e.g. gq — > gq from DIS where each flavor is weighted by electric 
charge squared. 

We choose A 1 as the object of the analysis, since it is more close to the direct observable 
in experiments than g^(x, Q 2 ). The g^(x, Q 2 ) data published by the experiments depend 
on the knowledge on the unpolarized structure functions at the time of their publication. 
By choosing A 1 as the object of the analysis, we can extend the analysis to include new 
set of data easily without any change in the previous data set. 

As explained in Sec. II, we parametrize the polarized parton distributions at small 
momentum transfer squared Q 2 = 1.0 GeV 2 (= Ql) with a special emphasis on the 
positivity and quark counting rule. Then, they are evolved to the Q 2 points, where the 
experimental data were taken, by the leading-order (LO) or NLO Q 2 evolution program. 
Using one of well-established unpolarized parton distributions, we construct A 1 as 



to compare with the experimental data. The polarized parton distributions at the initial 



PDFs. 




(1.1) 
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Ql are determined by a x 2 analysis. 

In Sec. II, we describe the outline of our analysis with the necessary formulation and 
the data set used in the analysis. Section III is devoted to the explanation of the LO and 
NLO Q 2 evolution programs which we developed for our fit. The parametrization of the 
polarized parton distribution functions at the initial Ql is described in Sec. IV, and the 
fitting results are discussed in Sec. V. The conclusions are given in Sec. VI. 



II. PARTON MODEL ANALYSIS OF POLARIZED DIS DATA 

In the experiments of polarized DIS, direct observables are the cross-section asymme- 
tries A\\ and A±, which are defined as 

A„ = ^ti, A ± = (2.1) 

The (7-f-f and a-^ represent the cross sections for the lepton-nucleon scattering with their 
parallel and anti-parallel helicity states, respectively. On the other hand, the cr-f_» and <7j__> 
are the scattering cross sections for transversely polarized nucleon target. We suppress 
the dependence on x and Q 2 where it is evident hereinafter. The asymmetries, A\\ and 
A±, are related to the photon absorption cross section asymmetries, Ai and A 2 , by 

A\\=V(A 1 + - n A 2 ), A ± = d(A 2 -(A 1 ), (2.2) 

where V represents the photon depolarization factor and rj is approximated as 7(1 — 
y)/(l — y/2) with 7 = 2Mx/ ' \[Q 2 . The d and ( are other kinematical factors. The 
asymmetries, A-y and A 2 , can be expressed as: 

A!(x,Q 2 ) 
A 2 (x,Q 2 ) 

Here a T 1 and a T 3 are the absorption cross sections of virtual transverse photon for the 

1 ' 2 > 2 

total helicity of the photon- nucleon system of \ and |, respectively; cjlt is the interference 
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<*r,i-<*r,§ g^x^Q^-^g^x^Q 2 ) 



(2.3) 



a T 1 + <j T 3 Fi (x, Q 2 ) 

' 2 '2 

"),-T... - [/;., r / Q 2 )] 



o" T i+o" T 3 Fi(x,Q 2 ) 



'2 '2 



term between the transverse and longitudinal photon-nucleon amplitudes; Fi(x, Q 2 ) is the 
unpolarized structure function of the nucleon. If we measure both An and A±, we can 
extract both gi(x, Q 2 ) and g 2 (x,Q 2 ) from experimental data with minimal assumptions. 
Otherwise, r/A 2 should be neglected in Eq. (2.2) to extract A 1 . This is justified since 
r)A 2 is much smaller than A\ in the present kinematical region. However, its effect has to 
be included in the systematic error. In the small-x or large-Q 2 region, r y 2 is the order of 
1CT 3 — 1CT 2 . An absolute value of g 2 (x, Q 2 ) has been measured to be significantly smaller 
than gi(x,Q 2 ). Therefore, the asymmetry in Eq. (2.3) can be expressed by 

to good approximation. Since the structure function usually extracted from unpolarized 
DIS experiments is F 2 (x,Q 2 ), we use F 2 (x, Q 2 ) instead of F ± (x, Q 2 ) by the relation 

F ^wv£m- (2 - 6) 

The function R(x, Q 2 ) represents the cross-section ratio for the longitudinally polarized 
photon to the transverse one, (Jl/ctt, which is determined experimentally in reasonably 
wide Q 2 and x ranges in the SLAC experiment of Ref. [18]. Recently published data 
on R(x, Q 2 ) by NMC [19] showed slightly different values from the SLAC measurement 
but mostly agreed within experimental uncertainties. Therefore, we decided to use SLAC 
measurements to be consistent with the most of the analyses of polarized DIS experiments. 

The structure function F 2 can be written in terms of unpolarized PDFs with coefficient 
functions as 

F 2 (x,Q 2 ) = ^e 2 xic q {x,a s ) <g> [ qi (x, Q 2 ) + q^x, Q 2 )] + C g (x, a,) <g> g(x, Q 2 ) \. (2.7) 
i=i ^ > 

Here qi and ^ are the distributions of quark and antiquark of flavor i with electric charge 
ej. The gluon distribution is represented by g(x,Q 2 ). The convolution ® is defined by 



f(x) ® 
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The coefficient functions, C q and C g , are written as a series in a s with x-dependent 
coefficients: 

CO , 

The LO coefficient functions are simply given by 

Cf\x) =8(l-x), C^(x)=0. (2.10) 
In the same way, the polarized structure function gi(x, Q 2 ) is expressed as 

9l (x, Q 2 ) = ^J2 e i\ AC l( X > "») [ A *( X ' Q ^ + A *( X ' Q ^ + AC ^ ^ A ^ X ' Q ^ f ' 
i=l ^ ' 

(2.11) 

where Ag^ = q\ — q\ {i = u, d, s, ...) represents the difference between the number densities 
of quark with helicity parallel to that of parent nucleon and with helicity anti-parallel. 
The definitions of Aqi and Ag are the same. The polarized coefficient functions AC q and 
AC g are defined similarly to the unpolarized case. 

Another separation of the quark distribution can be done by using flavor- singlet quark 
distribution AS(x, Q 2 ) and flavor- nonsinglet quark distributions for the proton and the 
neutron, Aq p NS (x, Q 2 ) and Aq r ^ s (x, Q 2 ), respectively. Those can be expressed with polar- 
ized PDFs as follows: 

AE(x) = a (x) = Au + (x) + Ad + (x) + As + (x), 
3 1 

= ±jj [ Au + (x) - Ad + (x) ] + i [Au + (x) + Ad + {x) - 2As + {x) ], (2.12) 

where Au + (x) = Au(x) + Au(x) and similarly for Ad + (x) and As + (x). Analyses in 
Ref. [12] and Ref. [15] utilized this separation. Such separation is useful in Q 2 evolution, 
and it is also natural when one wants to obtain quark contribution to the proton spin, 
Jo AE(x)dx. 

On the other hand, when we try to calculate the cross section for polarized pp reaction, 
e.g. Drell-Yan production of lepton pairs, we need the combination of Agj(xi) x Ag^a^) 
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(multiplied by electric charge squared). To allow such calculations with the above sepa- 
ration, we need further assumption on the polarized antiquark distributions, e.g. flavor 
symmetric sea, Au sea (x) = Au{x) = Ad sca (x) = Ad(x) = As(x) = As(x). With such 
assumption, the above separation becomes equivalent to the PDF separation in a sense 
that one description can be translated to another by simple transformation. 

Of course, we already know that unpolarized sea-quark distributions are not flavor 
symmetric [20] from various experiments including Drell-Yan production of lepton pairs in 
pp and pd collisions. Therefore, this assumption is only justified as an approximation due 
to limited experimental data. In principle, charged-hadron production data could clarify 
this issue. Although a x 2 analysis for the SMC and HERMES data seems to suggest a 
slight Am excess over Ad [21], the present data are not accurate enough for finding such 
a flavor asymmetric signature. Future experiments with charged current at RHIC [22] 
and polarized option at HERA will be very useful in improving our knowledge on the 
spin-flavor structure of the nucleon. Furthermore, as it has been done in the unpolarized 
studies, the difference between the polarized pp and pd cross sections provides a clue for 
the polarized flavor asymmetry [23] although actual experimental possibility is uncertain 
at this stage. 

The parametrization models studied so far have various differences in other aspects: 
(a) the choice of the renormalization scheme, (b) the functional form of the polarized 
parton distributions due to different physical requirements at Qq, and (c) the physical 
quantity to be fitted. In the following, we describe our position on these issues. 

• Renormalization Scheme 

Although the parton distributions have no scheme dependence in the LO, they do 
depend on the renormalization scheme in the NLO and beyond. In the polarized 
case, we have different choices of the scheme due to the axial anomaly and the 
ambiguity in treating the 75 in n dimensions [16]. In the NLO analysis, the widely- 
used scheme is the modified minimal subtraction (MS) scheme, in which the first 
moment of the nonsinglet distribution is Q 2 -independent. It was used, for example, 

8 



by Mertig and van Neerven [8] and Vogelsang [9]. However, the first moment of the 
singlet distribution is Q 2 -dependent in this scheme and thus it is rather difficult to 
compare the value of AT,(x,Q 2 ) extracted from the DIS at large Q 2 with the one 
from the static quark model at small Q 2 . To cure this difficulty, Ball, Forte and 
Ridolfi [24] used the so-called AB (Adler-Bardeen) scheme, in which the first moment 
of the singlet distribution becomes independent of Q 2 because of the Adler-Bardeen 
theorem [25]. In those schemes, however, some soft contributions are included in the 
Wilson coefficient functions and not completely absorbed into the PDFs. Another 
scheme called the JET scheme [26] or the CI (chirally invariant) scheme [27] has 
been recently proposed. All the hard effects are absorbed into the Wilson coefficient 
functions in this scheme. 

Although we choose the MS scheme in our analysis, the polarized PDFs in one 
scheme are related to those in other schemes with simple formulae [16]. 

• Functional Form of polarized PDF and Physical Requirements 

Different functional forms have been proposed so far for the polarized PDFs by 
taking account of various physical conditions. We choose the functional form with 
the special emphasis on the positivity condition and quark counting rule [28] at 



The positivity condition is originated in a probabilistic interpretation of the parton 
densities. The polarized PDFs should satisfy the condition 



This is valid in the LO since we can have the complete probabilistic interpretation 
for each polarized distribution only at the LO. Even in NLO, however, the positivity 
condition for the polarized cross section Aa with the unpolarized cross section a, 



Ql 



1.0 GeV 2 . 



Af l {x,Ql)\<f t {x,Ql). 



(2.13) 



Aa\<a, 



(2.14) 
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should still apply for any processes to be calculated with the polarized PDFs to the 
order of 0(a s ). Since it is very difficult to calculate the polarized and unpolarized 
cross sections of the NLO for all the possible processes, it is not realistic to determine 
the polarized NLO distributions by the positivity condition of Eq. (2.14). In our 
analysis, we simply require that Eq. (2.13) should be satisfied in the LO and also 
NLO at Q\. It is shown in Ref. [29] that the NLO Q 2 evolution should preserve the 
positivity maintained at initial Q^. 

In many cases, Regge behavior has been assumed for x — > 0, and the color coherence 
of gluon couplings has been also used at x ~ [30]. Furthermore, it is an interesting 
guiding principle that the polarized distributions have a similar behavior to the 
unpolarized ones in the large-x region [16]. Since behavior of the distributions at 
large x is determined by the term (1 — x) 13 in the functions, where (5 is a constant, 
we simply require that the polarized distributions should have the same (1 — x) 13 
term as the unpolarized ones. 

Those physical requirements and assumptions have to be tested by comparing with 
the existing experimental data. 

As for the choice of Qjj, it has to be large enough to apply perturbative QCD, but 
it should be small enough to maintain a large set of experimental data. We find 
Ql=1.0 GeV 2 to be a reasonable choice in our analysis. 

• Physical Quantities to be Fitted 

In most of the polarized experiments, the data have been presented for Ai(x,Q 2 ) 
and g\(x, Q 2 ). Some analyses [12,14,15,30] used the g±(x, Q 2 ) as data samples, while 
others [10,11,13,16] used the Ai(x, Q 2 ). It should be, however, noted that g\(x, Q 2 ) is 
obtained by multiplying A±(x, Q 2 ) by F±(x, Q 2 ), so that it is not free from ambiguity 
of the unpolarized structure function, Fi(x,Q 2 ). Therefore, we consider that it is 
more advantageous to use the A±(x, Q 2 ) as the data samples not only for the current 
work but also for the convenience in expanding the data set to include new data set 
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from SLAC, DESY (German Electron Synchrotron), CERN, and RHIC. 

Another important quantity which we should carefully consider is the cross sec- 
tion ratio R(x,Q 2 ) = (? L /cr T , where cr L and a T are absorption cross sections of 
longitudinal and transverse photons, respectively. In principle, nonzero R(x, Q 2 ) is 
originated from radiative corrections in perturbative QCD, higher twist effects, and 
target mass effects. Higher twist contribution to R(x, Q 2 ) is expected to be small 
in the large Q 2 region. So far, some analyses employed nonzero R(x,Q 2 ), while 
other analyses assumed R(x, Q 2 ) = 0. However, the latter is not consistent with the 
experimental analysis procedure, since R(x, Q 2 ) is also used for the evaluation of 
photon depolarization factor V. Indeed our analysis shows that world data prefer 
R(x, Q 2 ) 7^ 0: the x 2 increases significantly with R = 0. Therefore, we use nonzero 
R(x, Q 2 ) in fitting the data of A ± (x, Q 2 ). 

Table I summarizes experiments with published data on the polarized DIS [1-5,31-34]. 
These measurements cover a wide range of x and Q 2 with various beam species and 
energies and various types of polarized nucleon target (not shown in the table). The 
listed are the number of data points above Q 2 = 1.0 GeV 2 , and the total number of data 
points are 375. 

We use the data with minimal manipulation to analyze them in our framework so 
as to be consistent with the Q 2 evolution, the unpolarized parton distributions, and the 
function R(x, Q 2 ). For example, the E143 provides the proton data which are obtained by 
combining the results of different beam energies using the weights based on the unpolarized 
cross sections [4] (28 points), in addition to "rati/' data for each beam energy (81 points at 
Q 2 >1 GeV 2 ). Such weights depend on the choice of the unpolarized structure functions, 
which are being updated. To localize dependence on the unpolarized structure functions 
in the final manipulation for getting gi(x, Q 2 ), i.e. Ai(x, Q 2 ) multiplied by F±(x, Q 2 ), we 
decided to use the "raw/' data in our analysis. 

Table I also includes analysis methods. One of the major differences in the analysis 
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is the treatment of the A 2 (x,Q 2 ) and g 2 (x : Q 2 ) contributions to the gi(x : Q 2 )/Fi(x : Q 2 ). 
Some of SLAC experiments measured both Am and A±_ to enable direct extraction of g\jF\ 
and g 2 jF\. Other experiments included possible contribution of r\A 2 in their estimation 
of systematic errors. 

As mentioned above, the choice of the function R(x, Q 2 ) potentially affects Ai(x, Q 2 ), 
thus final results on polarized PDFs, since the function affects the photon depolarization 
factor V. While it was assumed to be constant in the analyses of the early days, its 
x-dependence and Q 2 -dependence have been found to be significant [18]. To reflect the 
most updated knowledge of R(x, Q 2 ) on our analysis, we have reevaluated the E130 and 
EMC data by using R 19go (x, Q 2 ) [18], which most of the experiments employed. However, 
we found changes of a few percent in EMC data and about 10% in E130 data: both of 
them are smaller than experimental errors. 

III. Q 2 EVOLUTION 

In our framework and in most of the analyses of structure functions in the parton 
model, the polarized parton distributions are provided at certain Q 2 {= Qq) with a num- 
ber of parameters, which are determined so as to fit polarized experimental data. The 
experimental data, in general, range over a wide Q 2 region. The polarized parton dis- 
tributions have to be evolved from Qq to the Q 2 points, where experimental data were 
obtained, in the x 2 analysis. In calculating the distribution variation from Qq to given Q 2 , 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations are used. 

To compare our parametrization with the data, we need to construct A ± (x, Q 2 ) from 
the polarized and unpolarized PDFs. Since the determination of the unpolarized PDFs 
is not in our main scope, we decided to employ one of the widely-used set of PDFs. 
Although there are slight variations among the unpolarized parametrizations, the calcu- 
lated F 2 (x,Q 2 ) structure functions are essentially the same because almost the same 
set of experimental data is used in the unpolarized analyses. The Gliick-Reya-Vogt 
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(GRV) unpolarized distributions [35] have been used in our analyses; however, the 
parametrization results do not change significantly even with other unpolarized distri- 
butions. We checked this point by comparing the GRV F 2 (x, Q 2 ) structure function with 
those of MRST (Martin-Roberts-Stirling-Thorne) [36] and CTEQ (Coordinated Theoreti- 
cal/Experimental Project on QCD Phenomenology and Tests of the Standard Model) [37] 
at Q 2 =5 GeV 2 in the x range 0.001 < x < 0.7. The differences between these distributions 
are merely less than about 3%. The differences depend on the x region; however, we find 
no significant systematic deviation from the GRV distribution. 

We calculate the GRV unpolarized distributions at Qq=1 GeV 2 in Ref. [35] 1 . The dis- 
tributions are evolved to those at Q 2 by the DGLAP equations, then they are convoluted 
with the coefficient functions by Eq. (2.7). Because the unpolarized evolution equations 
are essentially the same as the longitudinally polarized ones in the following, except for 
the splitting functions, we do not discuss them in this paper. The interested reader may 
read, for example, Ref. [38]. 

The polarized PDFs are provided at the initial Q$; therefore, they should be evolved 
to Q 2 by the DGLAP equation in order to obtain gi(x,Q 2 ). The DGLAP equations 
are coupled integrodifferential equations with complicated splitting functions in the NLO 
case. Both the LO and NLO cases can be handled by the same DGLAP equation form; 
however, the NLO effects are included in the running coupling constant a s (Q 2 ) and in 
the splitting functions APij(x). 

In solving the evolution equations, it is more convenient to use the variable t defined 

by 

t = \nQ 2 , (3.1) 



Actual calculation has been done by the FORTRAN program, which was obtained from the 
www site, http://durpdg.dur.ac.uk/HEPDATA/PDF. 
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instead of the variable Q 2 . Then, the flavor nonsinglet DGLAP equation is given by 

^Aq NS (x,t) = ^AP q±!NS (x)®Aq NS (x,t), (3.2) 

where Aq NS (x, t) is a longitudinally-polarized nonsinglet parton distribution, and AP q ± tNS 
is the polarized nonsinglet splitting function. The notation in the splitting function 
indicates a "Ag± Aq type" distribution ^ a^A^i Aqj), where Oj is given constant with 
flavor i. The singlet evolution is more complicated than the nonsinglet one due to gluon 
participation in the evolution. The singlet quark distribution is defined by AE(x,t) = 
Y^i ! {Aqi+ Aqi) , and its evolution is described by the coupled integrodifferential equations, 

-AK,„(.r) _\./-'„„(.r) ,A!Aj\ / I 

(3.3) 



d_ 

di 



AE(x,t) 
Ag(x,t) 



a s (t) 



2tt 



AP qq (x) AP qg (x) 
yAP gq (x) AP gg (x) ) yAg(x,t) ^ 

The numerical solution of these integrodifferential equations is obtained by a so-called 
brute-force method. The variables t and x are divided into small steps, 5U and 5xi 



respectively, and then the integration and differentiation are defined by 

df(x) f(x m +i) - f(x m ) 



dx 



Sx r 



, N x 

/ f(x) dx =^ Sx m f(Xr, 
J m=l 



(3.4) 
(3.5) 



The evolution equation can be solved numerically with these replacements in the DGLAP 
equations. This method seems to be too simple; however, it has an advantage over 
others not only in computing time but also in future applications. For example, the 
evolution equations with higher-twist effects cannot be solved by orthogonal polynomial 
methods. It is solved rather easily by the brute-force method [38]. Another popular 
method is to solve the equations in the moment space. However, the x distributions are 
first transformed into the corresponding moments. Then, the evolutions are numerically 
solved. Finally, the evolved moments are again transformed into the x distributions. If the 
distributions are simple enough to be handled analytically in the Mellin transformation, it 
is a useful method. However, if the distributions become complicated functions in future 
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or if they are given numerically, errors may accumulate in the numerical Mellin and inverse 
Mellin transformations. Therefore, our method is expected to provide potentially better 
numerical solution although it is very simple. 

The employed method is identical to that in Ref. [38] in its concept, but we had to 
improve the program in its computing time, since the evolution subroutine is called a 
few thousand times in searching for the optimum set of polarized distributions. There 
are two major modifications. The first one is to change the method of the convolution 
integrals, and the second is to introduce the cubic Spline interpolation for obtaining the 
parton distributions during the evolution calculation. Previously we calculated the con- 
volution integral by ^ AP(x/y)Aq(y, t). In this case, we had to calculate the splitting 
functions for each x value in the numerical integration, since the integration variable and 
the argument of the splitting function are different. Because the NLO splitting functions 
are complicated, this part of calculation consumed much time. In the present program, 
we evaluate the integral by ^ AP(y)Aq(x/y, t), which is mathematically equivalent to 
the above integral, and thus, we only need to calculate the splitting functions at a fixed 
set of x values once before the actual evolution. For example, the nonsinglet equation, 
Eq. (3.2), becomes 

Aq NS (x k ,t j+1 ) = Aq Ng (x k ,tj) + 8tj — AP q ± tNS (x m ) Aq NS ( — ,tA . (3.6) 

If the initial distribution Aq NS (xk/x m , to = 0) is provided, the next distribution 
Aq NS (xk, ti) is calculated by the above equation. Then, Aq NS (xk/x m ,ti) is calculated 
by the cubic Spline interpolation. Repeating this step N t — 1 times, we obtain the evolved 
nonsinglet distribution Aq NS (xk,tN t ). With these refinements, the evolution equations are 
solved significantly faster, and the subroutine can be used in the parametrization study. 

We show the Q 2 dependence in g^(x,Q 2 ) and A\(x, Q 2 ) as a demonstration of the 
performance of our program. The numerical calculations are done such that the accuracy 
becomes better than about 2% in the asymmetry A\. The LO and NLO (set NLO-1) 
parton distributions obtained in our analyses are used. The details of these distributions 
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are discussed in Sec. V. The initial structure functions g\ at Q 2 = 1.0 GeV 2 are evolved 
to those at Q 2 = 60.0 GeV 2 . Most of the used A\ data are within this Q 2 range. The LO 
and NLO results are shown in Fig. 1 by the dashed and solid curves, respectively. The 
LO distributions tend to be shifted to the smaller x region than the NLO ones. There 
are two reasons for the differences between the LO and NLO distributions. One is the 
difference between the LO and NLO F 2 structure functions for fitting the same data set 
of A 1 , and the other is the difference in Q 2 evolution. 

In Fig. 2, our Q 2 evolution curves at x — 0.117 are shown with the asymmetry A\ data 
by the SMC [2], SLAC-E143 [4], and HERMES [5] collaborations. The initial distributions 
are our LO and NLO parametrizations at Q 2 = 1 GeV 2 . The dashed and solid curves 
indicate the LO and NLO evolution results, respectively. In the large Q 2 region, both 
Q 2 variations (dAi/dlnQ 2 ) are almost the same; however, they differ significantly at 
small Q 2 , particularly in the region Q 2 < 2 GeV 2 . As the Q 2 becomes smaller, the 
NLO contributions become more apparent. We find that the theoretical asymmetry has 
Q 2 dependence although it is not large at x — 0.117. It is often assumed that the 
experimental asymmetry Ai(x, Q 2 ) is independent of Q 2 by neglecting the Q 2 evolution 
difference between gi(x, Q 2 ) and F ± (x, Q 2 ) in extracting the gi(x, Q 2 ) structure functions. 
The assumption has no physical basis. For a precise analysis, the Q 2 dependence in the 
asymmetry has to be taken into account properly and our framework is ready for such 
precision studies. 

IV. PARAMETRIZATION OF POLARIZED PARTON DISTRIBUTIONS 

Now, we explain how the polarized parton distributions are parametrized. The unpo- 
larized PDFs fi{x,Ql) and polarized PDFs Afi(x, Qq) are given at the initial scale Q^. 
Here, the subscript % represents quark flavors and gluon. These functions are generally 
assumed to be in a factorized form of a power of x inspired by Regge-like behavior at small 
x, a polynomial of x at medium x, and a power of (1 — x) expected from the counting 
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rule at large x: 

fi(x,Ql) = C iX a » (l-x) a ^ (1 +J2» 3i , j x a ^), (4.1a) 

j 

Af^Ql) = A jA' (1 - x)P" (1 • Yl ; >- ( 4 - lb ) 

where Cj and Dj are normalization factors and an, a^i, ot3i,j, cxuj, flu, fin, foij, and 
are free parameters. 

From the best fit to all the experimental data of the polarized DIS including new 
data, we can determine, in principle, the parameters in Eq. (4.1b). In practice, however, 
some of the parameters highly correlate each other and it is difficult to determine all the 
parameters independently. Therefore, it is desirable to reduce the number of parameters 
by applying physical conditions instead of leaving all these parameters free. 

In the present analysis, to constrain the explicit forms of polarized PDFs, we require 
two natural conditions: (i) the positivity condition of the PDFs and (ii) the counting rule 
for the helicity-dependent parton distribution functions. 

In order to make the positivity condition of Eq. (2.13) be tractable in the numerical 
analysis, we modify the functional form of the polarized PDF as 

*fi(x,Ql)=h i (x)f i (x,Q 2 ), (4.2) 

where 

hi{x) = A t x a * (1 -x) ft (1 + 7i x Al ), (4.3) 

at the initial scale Qq. Therefore, the positivity condition can be written as 

I hi(x) | < 1 (4.4) 

Furthermore, taking into account of the counting rule mentioned in Section II, we 
reduce Eq. (4.3) to 

h l {x)=A l x^{l+ li x x *), (4.5) 
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and we have the following functional form of polarized PDFs at Qfo: 

Afi(x,Ql) = Ai x ai (1 + 7 ^ A fi(x,Ql). (4.6) 

Thus, we have four parameters (Ai, ai, 7$ and A,) for each i. 

We further reduce the number of free parameters by assuming the SU(3) flavor sym- 
metry for the sea-quark distributions at Qq. As mentioned in Section II, this is simply 
a compromise due to a lack of experimental data. It should be noted that the sea-quark 
distributions are not SU(3) flavor symmetric at Q 2 > Ql even with the symmetric distri- 
butions at the initial Qq. 

When we assume this SU(3) flavor symmetric sea, the first moments of Au v (x) and 
Ad v (x) for the LO, which are written as rj Uv and i] dv , respectively, can be described in 
terms of axial charges for octet baryon, F and D measured in hyperon and neutron 
/3-decays as follows: 

Vu v ~Vd v = F + D, 

r) Uv +Vd v =3F-D. (4.7) 

Note that Eq. (4.7) is also used for the NLO (MS) case. Recently, since the /5-decay 
constants have been updated [39], we reevaluate F and D from the x 2 fit to the experi- 
mental data of four different semi-leptonic decays: n ^ p, A ^ p, S — > A, and £ — > n, by 
assuming the SU(3) / symmetry for the axial charges of octet baryon. With x 2 /d.o.f=0.98, 
the F and D are determined as 

F = 0.463 ±0.008, 

D = 0.804 ± 0.008, (4.8) 

which lead to rj Uv = 0.926 ±0.014 and r/ dv = —0.341 ±0.018. In this way, we fix these two 
moments at their central values, so that two parameters A Uv and A dv are determined by 
these first moments and other parameter values. Thus, the remaining job is to determine 
the values of remaining 14 parameters, Aq, A g , ai, 7$, Aj (i = u v , d v , q, g), by a x 2 
analysis of the polarized DIS experimental data. 
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V. NUMERICAL ANALYSIS 



A. x 2 analysis 

We determine the values of 14 parameters from the best fit to the Ai(x, Q 2 ) data for 
the proton (p), neutron (n) and deuteron (d). Using the GRV parametrization for the 
unpolarized PDFs at the LO and NLO [35] and the SLAC measurement of R(x,Q 2 ), we 
construct A c ^ lc (x, Q 2 ) for the p, n, and d. For the deuteron, we use gf = — \wd) 

with the D-state probability in the deuteron uj d = 0.05. 

Then, the best parametrization is obtained by minimizing x 2 — '^2(Af a,ta, (x, Q 2 ) ~ 
Al alc (x,Q 2 )) 2 /(AAf ata (x,Q 2 )) 2 with Minuit [40], where AAf ata represents the error on 
the experimental data including both systematic and statistical errors. Since some of 
the systematic errors are correlated, it leads to an overestimation of errors to include all 
systematic errors. On the other hand, if we fully exclude them, the uncertainties in the 
experimental data are not properly reflected in the analysis. Because of our choice to 
include the systematic errors, the x 2 defined in our analysis is not properly normalized. 
The minimum x 2 divided by a number of degree-of-freedom achieved in the analysis is 
often smaller than unity. Consequently the x 2 i n our analysis should be regarded as only 
a relative measure of the fit to the experimental data. In addition, the parameter errors 
are overly estimated. We have confirmed that inclusion of only statistical errors in the x 2 
analysis does not change the results significantly except a change of the x 2 by 7%, which 
is consistent with the change of the error size. 

In evolving the distribution functions with Q 2 , we neglect the charm-quark contribu- 
tions to Ai(x, Q 2 ) and take the flavor number Nf — 3 because the Q 2 values of the A\ 
experimental data are not so large compared with the charm threshold. To be consistent 
with the unpolarized, we use the same values as the GRV, A^ D = 204 MeV at LO and 
^■qcd = 299 MeV at NLO in the MS scheme. The NLO scale parameter leads to the value 
of a s (M§) = 0.118. In order to obtain a solution which satisfies the positivity condition, 
we make further refinements to the parametrization functions hi(x). The technical details 
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are discussed in Appendix A. 

The results are presented in Table II for the LO with x 2 /d.o.f=322. 6/360 and in Table 
III for the NLO with x 2 /d.o.f=300. 4/360. Because the values of A, are determined by the 
first moments for the Au v and Ad v distributions, they are listed without errors. We show 
the LO and NLO fitting results for the asymmetry A\ together with experimental data in 
Fig. 3. The theoretical curves are calculated at Q 2 =5 GeV 2 . The asymmetries are shown 
for the (a) proton, (b) neutron, and (c) deuteron. As the experimental data, the E130, 
E143, EMC, SMC, and HERMES proton data are shown in Fig. 3(a); the E142, E154, 
and HERMES neutron data are in (b); the E143, E155, and SMC deuteron data are in (c). 
Kinematical conditions and analysis methods of these experiments are listed in Table I. 
We find from these figures that the obtained parameters reproduce well the experimental 
data of A 1 in both LO and NLO cases. However, there are slight differences between the 
LO and NLO curves in Fig. 3, and three factors contribute to the differences. First, the 
most important difference is the contribution of the polarized gluon distribution through 
the coefficient function. Second, the LO and NLO evolutions are different because not 
only the the splitting functions but also the scale parameters are different. Third, the LO 
and NLO expressions are different in the unpolarized GRV distributions. 

B. Comparison of LO and NLO analyses 

Comparing the value of x 2 /d.o.f for the LO with that for the NLO, we found a better 
description of the experimental data with the NLO analysis. The value of x 2 /d.o.f is 
improved by 7%. This implies that it is necessary to analyze the data in the NLO if one 
wants to get better information on the spin structure of the nucleon from the polarized 
DIS data. 

The x 2 contribution from each data set is listed in Table IV. The improvement is 
significant especially for the HERMES proton and E154 neutron data. The results of g 1 
at the LO and NLO are shown in Fig. 4 and Fig. 5, respectively. The "experimental" 
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gi data are calculated by using Eqs. (2.5) and (2.6) together with the raw data for the 
asymmetry A\ and the GRV unpolarized distributions. The theoretical results are shown 
by the dashed, solid, and dotted curves at Q 2 =l, 5, 20 GeV 2 . As already shown in Fig. 
1, the #1 structure function shifts to the smaller-x region as Q 2 increases. It is rather 
difficult to discuss the agreement with the deuteron data in Figs. 4(c) and 5(c) because 
of the large experimental errors. However, the proton and neutron data at small x tend 
to agree with the theoretical curves at Q 2 =l GeV 2 . It is particularly clear in the neutron 
gi in Figs. 4(b) and 5(b). Furthermore, the proton, neutron, and deuteron data at large 
x agree with the LO and NLO curves at Q 2 =20 GeV 2 . There are correspondences of the 
data to the theoretical results because the small-x data are typically in the small-Q 2 range 
(Q 2 = 1 ~a few GeV 2 ) and the large-x data are in the large-Q 2 range (Q 2 >10 GeV 2 ). 

As seen in Figs. 4 and 5, the LO gf is slightly larger at small x in comparison with the 
NLO gf, while the LO g? is smaller than the NLO g? in the range 0.01 < x < 0.2. The 
NLO fit agrees better with the data. The x 2 improvement in the NLO for the HERMES 
and E154 data in Table IV is explained as follows by using Fig. 3. In comparing the 
theoretical curves with the data, we should note that the theoretical asymmetries are 
given at fixed Q 2 (Q 2 =5 GeV 2 ), whereas the data are at various Q 2 values. However, as it 
is found in Fig. 3(a), the LO curve is slightly above the NLO one and also the HERMES 
data. It makes the x 2 value larger in the LO analysis. In Fig. 3(b), it is clear that 
the LO curve deviates from the E154 neutron data, so that the x 2 contribution becomes 
larger from the E154 data. It is well known that the difference between the NLO (MS 
scheme) and LO originates from the polarized gluon contribution to the structure function 
g 1 via the Wilson coefficient. Accordingly, the result that the NLO fit is better than the 
LO implies that the polarized gluon has a nonzero contribution to the nucleon spin, i.e. 
Ag at Ql- Furthermore, we find in this analysis that the NLO fit is more sensitive to 
the polarized gluon distribution than the LO one. Therefore, we can conclude that the 
NLO analysis is necessary to extract information on the polarized gluon distribution. 
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C. Behavior of polarized parton distribution functions 

We show the behavior of polarized parton distributions xAfi(x,Q 2 ) as a function of 
x at Q 2 = 1 GeV 2 for the (a) LO and (b) NLO cases in Fig. 6. The first moment for 
Au v (x) is fixed at the positive value (i] Uv =0.926) and the one for Ad v (x) is at the negative 
value (i] dv = —0.341), so that the obtained distributions Au v (x) and Ad v (x) become 
positive and negative, respectively. In the same way as the other x 2 -analysis results, 
the antiquark (gluon) distribution becomes negative (positive) at small- and medium- x 
regions. The gluon distribution cannot be determined well by only the lepton scattering 
data. In particular, the gluon distribution plays a role in g\ only through the Q 2 evolution 
in the LO, so that Ag(x) cannot be uniquely determined. Even if it is neglected in the 
analysis (Ag = 0), the x 2 difference is not so significant in the LO. The NLO effects are 
apparent by comparing Fig. 6(a) with Fig. 6(b). In the NLO, the gluon distribution 
contributes to g\ additionally through the coefficient function; therefore, it modifies the 
valence-quark distributions (particularly the Au v ) and the antiquark distribution. The 
NLO distribution Au v becomes significantly smaller than the LO one at small x, and 
the NLO distribution Aq becomes a more singular function as x — > 0. Because of more 
involvement of the gluon distribution in gi, the determination of Ag is better in the NLO 
X 2 analysis. 

Recently, the measurement of polarized parton distributions of each flavor has been 
carried out by the SMC in semi-inclusive processes of the polarized DIS [41]. Although we 
did not include the semi-inclusive date in our analysis from the consideration of the data 
precision and the analysis framework, it is still possible to compare our polarized PDFs 
with their analysis. In order to compare with the SMC data, the LO initial distributions 
are evolved to those at Q 2 = 10 GeV 2 by the LO evolution equations. Then, the ratios 
Au v (x)/u v (x) and Ad v (x)/d v (x) are shown in Fig. 7 together with the SMC data. The 
theoretical ratios are roughly constants in the small-x region (x < 0.1) and Au v (x)/u v (x) 
approaches +1 as x — > 1 whereas Ad v (x)/d v (x) approaches —1. We find that our LO 
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parametrization seems to be consistent with the data. However, it is unfortunate that 
our NLO parametrization cannot be compared with the data since the SMC data are 
analyzed only for the LO. 

D. Small- a; behavior of polarized antiquark distributions 

As we obtained in the \ 2 analyses, the small-x behavior of the parton distributions 
is controlled by the parameter a. It is obvious from Tables II and III that the small- a; 
behavior cannot be determined in the antiquark and gluon distributions. For example, 
the obtained parameter is listed as Oiq(NLO) = 0.32 ± 0.22 with a large error. It suggests 
that the small- a; part of the antiquark distribution cannot be fixed by the existing data. 
In order to clarify the situation, we need to have higher-energy facilities such as polarized- 
HERA and eRHIC [42]. 

Because the present experimental data are not enough for determining the small- a; 
behavior, we should consider to fix the parameter a for the antiquark distribution by 
theoretical ideas. The gluon parameter a g cannot be also determined. However, we leave 
the problem for future studies because the lepton scattering data are not sufficient for 
determining the gluon distribution in any case. Some predications are made for a.q in the 
following by using the Regge theory and the perturbative QCD. 

According to the Regge model, the structure function gi in the small- a; limit is con- 
trolled by the intercepts (a) of ai(1260), /i (1285), and /i(1420) trajectories: 

gi(x) ~ x~ a as x — > 0. (5.1) 

However, not only the a\ intercept but also the f\ intercepts are not well known. It is usu- 
ally assumed as a ai = —0.5 ~ [43]. Therefore, we expect Ag ~ 3.(0.0,0.5)^ w h ere 3.(0.0,0.5) 
indicates that the function is in the range from a; to x 5 . Since our parametrization is 
provided for the function hi(x) = Afi(x)/ fi(x), we should find out the small-x behavior of 
the unpolarized distribution. According to our numerical analysis, the GRV distribution 
has the property ig ~ x~°- 1A at Q 2 =l GeV 2 . Taking these small- a; functions into account, 
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the Regge prediction is 



hf" e {x)^x { - 1A ^\ 



(5.2) 



if the theory is applied at Q 2 = 1 GeV 2 . Our LO and NLO fits result in a; - 59 and x a32 , 
respectively, as x — > 0. These functions look very different from Eq. (5.2); however, they 
are not inconsistent if the errors of Tables II and III are taken into account. 

The perturbative QCD could also suggest the small-x behavior. In the small-x limit, 
the splitting functions are dominated by the most singular terms. Therefore, if we can 
assume that the singlet-quark and gluon distributions are constants at certain Q 2 (= Q 2 ) 
in the limit x — > 0, their singular behavior is predicted from the evolution equations. 
According to its results, the singlet distribution behaves like [44] 



where £(Q 2 ) = \n[a s (Qi) /a s (Q 2 )}, C A = 3, and (3 = 11 - 2N f /3. The problem is to 
find an appropriate Q\ where the singlet and gluon distributions are flat at small x. 
Choosing the range Q\ = 0.3 ~ 0.5 GeV 2 and Q 2 = 1 GeV 2 , we fit the above equation 
numerically by the functional form of x~ a at small x. Then, the obtained function is in 
the range, x^~ 0A2 ' -°- 09 ). Because the unpolarized distribution is given by xq ~ x~°- u , the 
perturbative QCD (with the assumption of the above Q\ range) suggests 



This function falls off much faster than ours at small x. 

In this way, we found that the perturbative QCD and the Regge theory suggest the 
small-x distribution as hq ~ ^(i- . 1 - 6 ). Because the small x behavior cannot be determined 
by the x 2 analyses in Sec. VA, we had better fix the power of x by these theoretical 
implications. In this subsection, the NLO \ 2 analyses are reported by fixing the parameter 
at ag=0.5, 1.0, and 1.6. The middle value is the perturbative QCD estimate, and the 
latter two ones are roughly in the Regge prediction range. The first one is taken simply 
by considering a slightly singular distribution than these theoretical predictions. 




(5.3) 




1.0 



(5.4) 
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The obtained parameters and x 2 are listed in Table V. Considering the NLO value 
x 2 =300.4 in Table IV, we find that the X 2 change is 0.1%, 1.8%, and 7.7% for a,=0.5, 1.0, 
and 1.6, respectively. The x 2 changes are so small in 0:^=0.5 and 1.0 that they could be 
equally taken as good parametrizations in our studies. Using the obtained distributions 
with fixed ag, we have the first moments and spin contents in Table VI. Because of 
the small- x falloff for larger etq, the antiquark first moment and spin content change 
significantly. If the perturbative QCD and Regge prediction range (ag=1.0 and 1.6) is 
taken, the calculated spin content is within the usually quoted values AS = 0.1 ~ 0.3. 
The obtained x 2 value suggests that the ag=1.0 solution could be also taken as one of 
the good fits to the data. In this sense, our results are not inconsistent with the previous 
analyses. However, the results indicate that a better solution could be obtained for smaller 
a^, so that the spin content could be smaller than the usual values AE = 0.1 ~ 0.3. At 
least, we can state that the present data are not taken at small enough x, so that the spin 
content cannot be determined uniquely. 

We found that the a^=0.5 and 1.0 results could be also considered as good parametriza- 
tions to the experimental data. The x 2 is so large in the ag—1.6 analysis that its set cannot 
be considered a good fit to the data. Because the 0^=0.5 results are almost the same as 
the NLO ones in Sec. V B, it is redundant to take it as one of our parametrizations. There- 
fore, we propose the LO and NLO distributions (sets: LO and NLO-1) in VB together 
with the ag=1.0 distributions (set: NLO-2) as three sets of the AAC parametrizations. 

Although the parametrization for A/j//j is necessary for imposing the positivity con- 
dition, it is rather cumbersome for practical applications in calculating other cross sections 
in the sense that we always need both our parametrization results and the GRV unpolar- 
ized distributions at Q 2 =l GeV 2 . Furthermore, it is not convenient that the analytical 
GRV distributions are not given at Q 2 =l GeV 2 . In Appendix B, we supply simple func- 
tions for the three AAC distributions without resorting to the GRV parametrization for 
the practical calculations. 



25 



E. Spin contents of polarized quarks and gluons 

The first moment of each polarized parton distribution and the integrated g\ at Q 2 =l, 
5, and 10 GeV 2 are given in Table VII for the LO and NLO. At Q 2 = 1 GeV 2 , the amounts 
of quarks and gluons carrying the nucleon spin are 

AS = 0.201, Ag = 0.831, in the LO, 

AS = 0.051, Ag = 0.532, in the NLO-1, (5.5) 
AS = 0.241, Ag = 0.533, in the NLO-2. 



These results confirm that the quarks carry a small amount of the nucleon spin. The first 
moments of the structure functions at Q 2 = 1 GeV 2 are 



r?(Q 2 ) 


= 0.144, 


I7(Q 2 ) = 


-0.067, 


TfiQ 2 ) 


= 0.036, 


in the LO, 




r?(Q 2 ) 


= 0.110, 


r7(Q 2 ) = 


-0.069, 


Tl{Q 2 ) 


= 0.019, 


in the NLO-1, 


(5.6) 


r?(Q 2 ) 


= 0.128, 


r7(Q 2 ) = 


-0.051, 


rf(Q 2 ) 


= 0.035, 


in the NLO-2. 





Because the first moment of Au v — Ad v is fixed by Eq. (4.7), the Bjorken sum rule is 
satisfied in both LO and NLO at any Q 2 within the perturbative QCD range. 

It should be noted that our AS in the NLO-1 seems to be considerably smaller than 
the usual values published so far in many other papers. In fact, the recent SMC and 
Leader-Sidrov-Stamenov (LSS) parametrizations [15,16] obtained AS =0.19 and 0.28, 
respectively, at Q 2 =l GeV 2 . The difference originates mainly from the small-x behavior 
of the antiquark distribution. We compared our NLO-1 distribution Ag, which is denoted 
as AAC, with the other MS distributions in Fig. 8. The LSS (1999) antiquark distribution 
is directly given in their parametrization, whereas the SMC distribution is calculated by 
using their singlet and nonsinglet distributions. Because the antiquark distribution is 
not directly given in the SMC analysis, we may call it as a transformed SMC ("SMC") 
distribution. The transformed SMC has peculiar x dependence at medium and large x; 
however, all the distributions agree in principle in the region (0.01 < x < 0.1) where 
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accurate experimental data exist and the antiquark distribution plays an important role. 
On the other hand, it is clear that our distribution does not fall off rapidly as x — > in 
comparison with the others. This is the reason why our NLO-1 spin content is significantly 
smaller. 

In order to clarify the difference, we plot the spin content in the region between x m i n 
and f by calculating AS(x miri ) = f 1 AT l (x)dx in Fig. 9. Because the LSS and SMC 
distributions are less singular functions of x, their spin contents saturate even at x — ICT 4 
although our AE still decreases in this region. The difference simply reflects the fact that 
the accurate experimental data are not available at small x. The parametrization results 
with fixed ctq are also shown. As the antiquark distribution becomes less singular, the 
spin content becomes larger. As mentioned in the previous subsection, the ag—LO results 
could be taken as a good fit. The spin content is 0.24 in this case and it is completely 
within the usual range AS = 0.1 ~ 0.3. 

The small- x issue has been discussed in other publications. The idea itself stems from 
the publication of Close and Roberts [45], and it is also noted in the numerical analyses 
of Altarelli, Ball, Forte, and Ridolfi (ABFR) [12]. In the ABFR parametrization, various 
fits are tried by assuming the small-a; behavior, and they obtain the first moment of 
ao(x) as ao = 0.02 ~ 0.18. Therefore, our NLO-1 analysis is consistent with their studies 
although the spin content seems to be smaller than the usual one (0.1 ~ 0.3). In this 
way, our NLO-1 analysis result may seem very different from many other publications, 
it is essentially consistent with them. It indicates that the small-x (~ 10~ 5 ) data are 
absolutely necessary for the determination of the spin content. 



F. Comparison with recent parametrizations 

We have already partially discussed the comparison with the recent parametrization 
results in the previous subsection. However, the detailed discussions are necessary par- 
ticularly on the differences between these analyses in order to clarify the difference in the 
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physical basis. 

First, we discuss differences between our parametrization and the LSS. Before the 
detailed comparison, we used their x 2 -fitting procedure in our program and confirmed 
their numerical results. It indicates that both fitting programs are consistent although 
evolution methods and other subroutines are completely different. 

Our parametrization functions are similar to theirs. In fact, both methods use 
the parametrization for the ratio of the polarized distribution to the unpolarized one 
(Afi(x)/fi(x) = hi(x), i = u v , d v , q, g). The LSS parametrization employed a very sim- 
ple function hi(x) = AiX a % and we used a more complicated one hi(x) = AiX ai (l +jiX Xi ). 
This may seem to be insignificant; however, the extra parameters provide wide room for 
the functions to readjust in the x 2 analysis. According to our studies, the minimum x 2 
cannot reach anywhere close to our minimum point if the LSS function is used in our fit. 
Therefore, although it is a slight modification, the outcome has a significant difference. 
Furthermore, the LSS gluon distribution fails to satisfy the positivity condition at large 
x although it does not matter practically at this stage. 

Another important difference is how to calculate the spin asymmetry A 1 from the 
unpolarized distributions. There are two issues in this calculation procedure. One is that 
LSS kept the factor 1 + AM^x 2 /Q 2 in handling the SLAC data, whereas we neglected. 
Another is that LSS calculated the structure function F± directly from the unpolarized 
distributions, whereas we calculated it by Eq. (2.6). As for the first point, we have 
checked that inclusion of the factor has not significant impact on the results. It is partly 
because the factor 1 + AM^x 2 /Q 2 modifies the asymmetry A 1 at large x but the Q 2 values 
are generally large in such a x region. The second point is more serious. Their method is 
right in the light of perturbative QCD. However, the F 2 structure functions are generally 
used rather than Fi in obtaining the unpolarized PDFs. If there were no higher-twist 
contributions, it does not matter whether F\ is calculated directly or Eq. (2.6) is used. 
However, it is well known that the higher-twist effects are rather large as obvious from the 
function R(x,Q 2 ) in the SLAC-1990 analysis [18]. It modifies the asymmetries as large 
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as 35%, and the modification is conspicuous in the whole x region. In the LSS analysis, 
perturbative QCD contributions to the function R are included due to the coefficient- 
function difference between F\ and F 2 , but they are small in the small- and medium- x 
regions. This difference in handling F 1 creates the discrepancy between the LSS and 
our polarized antiquark distributions, and it is especially important for determining their 
small-a; behavior. 

Next, we discuss comparison with the SMC parametrization. Our x 2 analysis is dif- 
ferent from theirs in the parametrization functions. We parametrized the ratios A/j//, 
(i = u v , d v , q, g). As mentioned in Sec. II, the analysis by the SMC in Ref. [15] utilized 
the separation of the polarized quark distributions into (AE(x), Aq^ s (x), and Aq$ s (x)) 
which can, in principle, be transformed into Au + (x), Ad + (x), and As + (x). 

When we do this transformation of SMC results to compare with the polarized sea- 
quark distributions from our analysis and LSS, we find that the polarized strange-quark 
distribution (As(x)) from the "transformed SMC" oscillates as shown in Fig. 8. However, 
this simply implies that the conventionally used functional form has a limitation and the 
distribution functions obtained from different separations can be quite different. The 
uncertainty of the sea-quark distribution was also pointed out in the analysis by Gordon, 
Goshtasbpour, and Ramsey [14]. We should re-emphasize that direct measurement of the 
sea-quark polarization is very important. At the highest energy of polarized pp collisions 
at RHIC, the weak bosons are copiously produced and the parity violating asymmetry 
Al for its production is very useful in elucidating spin-flavor structure of the nucleon [46]. 
With such direct measurement, the uncertainty in the polarized sea-quark distribution 
will be much reduced. 

Common differences from the SMC and LSS are that a large set of data tables is used 
for Ai rather than the Q 2 averaged one. Although the present data may not have the 
accuracy to discuss the Q 2 dependence, it is desirable to use the large table if one wishes 
to obtain better information on the gluon distribution. Furthermore, an advantage of our 
results is that the positivity condition is strictly satisfied, so that our parametrizations 
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does not pose any serious problem in practical applications. 



VI. CONCLUSIONS 

We have analyzed the experimental data for the spin asymmetry A 1 of the proton, 
neutron, and deuteron by using a simple parametrization for the ratios of polarized parton 
distributions to the corresponding unpolarized ones. We discussed the details on physical 
meanings behind our parametrization and also on our Q 2 evolution method. As a con- 
sequence, we found that the asymmetry A 1 could have significant Q 2 dependence in the 
small Q 2 region (Q 2 < 2 GeV 2 ), so that frequently-used assumption of the Q 2 indepen- 
dence in Ai cannot be justified in a precise analysis. From the LO and NLO x 2 analyses, 
we obtained good fits to the experimental data. Because the NLO x 2 is significantly 
smaller than that of LO, the NLO analysis should be necessarily used in the parametriza- 
tion studies. An advantage of our analysis is that the positivity condition is satisfied in 
the whole x region. An important consequence of our analyses is that the small-a; behav- 
ior of the sea-quark distributions cannot be uniquely determined by the present data, so 
that the usual spin content AS = 0.1 ~ 0.3 could be significantly modified depending on 
the future experimental data at small x (~ 10~ 5 ). Our LO and NLO analyses suggested 
AE=0.20 and 0.05, respectively. However, if we take theoretical suggestions by "pertur- 
bative QCD" and Regge theory for the polarized antiquark distribution at small x, the 
spin content becomes AS = 0.24 ~ 0.28 in the NLO. The obtained gluon distributions are 
positive in both LO and NLO, but it is particularly difficult to determine Ag in the LO. 
From these analyses, we have proposed one LO set and two NLO sets of parametrizations 
as the AAC polarized parton distribution functions. 
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APPENDIX A: TREATMENT OF POSITIVITY CONDITION 

IN OUR x 2 ANALYSIS 

Additional modification of the function hi(x) is desirable in the actual x 2 fitting. 
Although Eq. (4.5) is a useful functional form, it is not very convenient for the x 2 
analysis in the sense that the positivity condition is rather difficult to be satisfied. In 
fact, running our x 2 program, we obtain a solution which does not necessarily meet the 
positivity requirement. In order to take into account this condition, the function is slightly 
modified although it is equivalent in principle: 



where 5i = & + K{. It can be seen why this function is more suitable at x = 1 by 
the following simple example. The original function is given by two parameters, hi(x = 
1) = Ai (1 + 7i); however, the modified one is by only one parameter hi(x = 1) = S^. 
Therefore, it is more easier to restrict the function hi(x) within the positivity-condition 
range. There is another advantage that the parameters are rather independent each other. 
For example, the parameter Aj is strongly correlated with (Aj > — ctj) if we would like 
to avoid singular behavior as x — > 0. In this way, the functional form of Eq. (Al) is used 
in the actual x 2 fitting although it is mathematically equivalent to Eq. (4.5). 

Although we could perform the x 2 analysis with the supplied information, it is not 
straight forward to obtain a solution which satisfies the positivity condition. We de- 
scribe the details of the analysis procedure. First, it was already mentioned that the first 
moments of Au v and Ad v are fixed by the F and D values, and they are given by 



(Al) 



= 5iX n - Ki(x Vi - X^) 



i = u v , d v , q, g 




(A2) 
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Then, the parameters k Uv and Ka v are determined by 

Sj j dx x Vi fjjx) - rn 

Ki jdx(x^-x^)fi(xy 1 j 

As we explained in Sec. VD, theory suggests the functions hi should not be a singular 
function of x in the small-x region. Therefore, we try to find a solution in the parameter 
range fj, h v % > 0. 

Next, we discuss the positivity condition. If the signs of the parameters and Ki are 
the same, the function hi(x) is a monotonically increasing or decreasing function, so that 
hi(x = 1) = 5i should be within the range — 1 < 5i < +1 due to the positivity requirement. 
On the other hand, if the signs are different, the function could have an extreme value 
at certain x (= X). If X is larger than one, the function could be a monotonic one in 
the range (0 < x < 1). Then, the same condition — 1 < Si < +1 is applied. However, if 
X is smaller than one, the situation is slightly complicated. Because the first and second 
terms have the same functional form in the first equation of Eq. (Al), we can have either 
/ij < Vi or Hi> Vi. Therefore, the condition < Vi is taken (practically only for Aq and 
Ag) in the following analysis without loosing generality. From Eq. (Al), we find that the 
extreme value is located at 

where d = Hi/vi (0 < Q < 1). It is in the range < X < 1 if the condition < —KiCi/d < 
1, namely 

t^V < K4 for < K4, 

Wi (A5) 
> Ki for Ki < 0, 



i-C< 

is satisfied. The extreme value is then obtained as 

— 

hi(X) = ( _J^iCi j ^ ^ (1 _ (A6) 
V Oi- KiJ 

Using the positivity condition < 1, we obtain the following constraint on the 

parameters: 

g + (Ki) = Ki-5i- K^i " Q)] 1 ^ > 0, (A7) 
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in the case «j > (0 < hi(X) < 1). Because the function g + (Ki) has a positive curvature, 
we try to find a «j point (= k;^) which satisfies = 0. There is only one solution for 

negative 5i and two solutions for positive <5j. In any case, we seek the solution k\ which 
is larger than the extreme point Ki = 1/(1 — d) by the Newton's method. Then, the 
parameter Ki is redefined The parameters Oi are used in the x 2 analysis for 

the antiquark and gluon distributions within the range < <7j < 1, so that the actual 
functional form is 

hi(x) = 5iX ai - (TiK'i(x ai - x 01 ^ 1 ) for % = q, g. (A8) 

On the other hand, we find 

g-(Ki) = K t -5 t - k^ - d)]^ < 0, (A9) 

in the case k { < (—1 < hi(X) < 0). A similar analysis is done for the function <7~(/Cj) 
in order to satisfy the positivity condition. With these preparations, we can perform the 
X 2 analysis. 
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APPENDIX B: PRACTICAL POLARIZED PARTON DISTRIBUTIONS 



Our polarized parton distributions are given in the parametrized functions hi(x) mul- 
tiplied by the GRV unpolarized distributions. For practical applications, we supply the 
following three sets of simple functions, which reproduce the x 2 analysis results in Sec. 
V, as the A AC distributions at Q 2 =l GeV 2 : 

Set : AAC - LO 

xAu v (x) = 0.4949 x a456 (l - xf m (l + 9.60 x 1 - 23 ), 
xAd v {x) = -0.2040 :r a456 (l - xf- 77 (I + 14.6 :r L36 ), 
xAq(x) = -0.1146 rr a536 (l - x) ia5 (l + 39.4 x im ), 
xAg{x) = 2.738 x a908 (l - xf M (l + 12.3 a; 1 ' 60 ), (Bl) 

Set : AAC - NLO - 1 

xAu v (x) = 0.4029 x 0A78 (l -x) 3 - 18 (l + 15.1 x im ), 
xAd v (x) = -0.2221 x a568 (l -x) 3 - 92 (l + 9.46 x - 813 ), 
xAq(x) = -0.03249 x a230 (l - x) 7 - 77 (l + 3.65 x a883 ), 
xAg(x) = 8.844x L77 (l - xf- 2 \l + 13.6x L51 ), (B2) 

Set : AAC - NLO - 2 

xAu v {x) = 0.4353 x a465 (l - x) 2M (l + 8.98 x a938 ), 
xAd v (x) = -0.1850 x 0A71 {l - x) 3 - S9 {l + 14.0 a; 1 ' 11 ), 
xAq(x) = -0.2452 x a752 (l - x) 8 ' 13 , 

xAg{x) = 8.895 x 1 ' 77 (I - xf 22 {l + 13.6 x 1 ' 51 ). (B3) 
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TABLES 

TABLE I. Summary of published polarized DIS experimental data on the spin asymmetry 
A\. The listed are the number of data points above Q 2 =l GeV 2 . 



Exp. 


x range 


Q 2 range 


#ofdata A 2 /g 2 (x,Q 2 ) 


R(x,Q 2 ) 


Ref. 












contribution 






EMC (p) 


0.015- 


-0.466 


3.5-29.5 GeV 2 


10 


neglected 


(5 2 -indep 


L J 


SMC (p) 


0.005- 


-0.480 


0.25-72.07 GeV 2 


12 


neglected 


±(a oon 


[21 

L J 


E130 (v) 


0.18- 


-0.70 


3 5-10 GeV 2 


8 


neglected 


constant 


[311 


E143 (p) 


0.022- 


-0.847 


0.28-9.53 GeV 2 


81 


measured 


-TL1 QQO 


L J 


HERMES (p) 


0.021- 


-0.85 


0.8-20.0 GeV 2 


19 


E143/SMC 




[51 

L J 


SMC (d) 


0.005- 


-0.480 


1.3-54.4 GeV 2 


12 


neglected 




[21 

L J 


E143 (d) 


0.022- 


-0.847 


0.28-9.53 GeV 2 


81 


measured 


-Rl990 


[4] 


E155 (d) 


0.01- 


-0.9 


1.0-40.0 GeV 2 


24 




-Rl990/NMC 


[34] 


E142 (n) 


0.035- 


-0.466 


1.1-5.5 GeV 2 


8 


neglected 


-Rl990 


[32] 


E154 (n) 


0.0174- 


-0.5643 


1.21-15.0 GeV 2 


11 


measured 


-Rl990 


[33] 


HERMES (n) 


0.033- 


-0.464 


1.22-5.25 GeV 2 


9 


neglected 


-Rl990 


[5] 



TABLE II. Obtained parameters at Q 2 = 1 GeV 2 in the leading-order \ 2 analysis. 



distribution 


A 


a 


7 


A 


AUy 


0.404 ± 0.054 


0.00 ± 0.01 


1.47 ± 0.20 


1.41 ± 0.46 


Ad v 


-0.274 ± 0.056 


0.00 ± 0.01 


2.65 ± 0.54 


1.25 ± 0.28 


Aq 


-0.680 ± 0.373 


0.59 ± 0.94 


-2.47 ± 0.82 


4.06 ± 1.14 


&g 


47.5 ± 4.1 


1.44 ± 0.73 


-0.986 ± 0.002 


0.06 ± 1.05 
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TABLE III. Obtained parameters at Q 2 = 1 GeV 2 in the next-to-leading-order x 2 analysis 
(set NLO-1). 



distribution 


A 

A 


Ot 


7 


\ 

A 


Au v 


0.356 ± 0.047 


0.00 ± 0.00 


1.54 ± 0.20 


0.889 ± 0.058 


Ad v 


-0.502 ± 0.031 


0.153 ± 0.065 


0.992 ± 0.062 


2.48 ± 0.27 


Aq 


-0.269 ± 0.107 


0.32 ± 0.22 


-4.72 ± 1.48 


3.20 ± 0.47 


A 5 


249.2 ± 8.3 


2.15 ± 0.11 


-1.0040 ± 0.0002 


0.031 ± 0.152 



TABLE IV. x 2 contribution of experimental data compared with the number of data points. 
Here, the NLO indicates the set NLO-1 . 



experimental 


# of data 




x 2 




data 


(Q 2 >1.0 GeV 2 ) 


LO 




NLO 


EMC (p) 


10 


5.2 




4.6 


SMC (p) 


59 


55.0 




53.7 


E130 (p) 


8 


5.1 




5.2 


E143 (p) 


81 


65.0 




60.8 


HERMES (p) 


19 


23.1 




17.2 


SMC (d) 


65 


56.6 




54.0 


E143 (d) 


81 


79.1 




81.2 


E155 (d) 


24 


20.0 




17.1 


E142 (n) 


8 


3.5 




2.4 


E154 (n) 


11 


7.5 




1.8 


HERMES (n) 


9 


2.6 




2.3 


total 


375 


322.6 




300.4 
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TABLE V. Obtained parameters at Q 2 = 1 GeV 2 in the NLO \ 2 analysis with fixed aq. 



distribution 


2 

X 


A 

A 


a 


7 


A 


(dig = 0.5) 


300.7 










Au v 




0.357 ± 0.093 


0.000 ± 0.001 


1.55 ± 0.40 


0.900 zb 0.335 






— u.oiz it u.uyo 


n 1 kg j- n 99V 


n nc;9 i n i oi 
u.yoz ± u.ioi 


Z.DO zt U.DD 


A — 

Aq 




n cr n i i ;i m " 

—9.50 ± 10.07 


U.5 (fixed) 


n non i n / \ / • / \ 

—0.980 ± 0.060 


n m no i n non A 

U.0102 ± 0.0394 


A 




1 /( O r I w A 

148.5 ± 7.4 


oil i n o i • 

2.11 ± 0.26 


— 1.0067 zb 0.0005 


0.051 zb 0.337 


/ i n\ 

(a, = 1.0) 


305.8 










Au v 




0.589 ± 0.055 


0.120 zb 0.090 


0.632 zb 0.059 


1.62 zb 0.27 






— U.Z /y ± U.Uoo 


U.UUU zt U.UU1 


O £C _L n on 


i qo i n oi 
l.oz dz U.ol 


A — 

Aq 




—47.7 ± 9.7 


1.0 (nxed) 


— 1.0065 ± 0.0056 


U.02U4 ± U.0707 


A 




173.8 ± 17. 3 


2.14 zb 0.19 


— 1.0058 zb 0.0007 


0.045 zb 0.253 


(a, = 1.6) 


323.5 










La U v 




i.oou in u.ioz 


U.OOO zt U.1ZU 


— / zt U.U41D 


A Q1 Q ^ A 90Q 

u.oio zt u.zuy 


Ad v 




-0.321 ± 0.097 


0.000 ± 0.000 


2.12 zb 0.64 


2.21 zb 0.17 


Aq 




-119.6 ± 11.8 


1.6 (fixed) 


-0.9976 zb 0.0071 


A AQ 9 _1_ 1 QA/1 

U.U33 ± 1.3U4 






176.6 ± 9.3 


2.77 ± 0.38 


-1.0057 zb 0.0004 


0.057 zb 0.547 


TABLE VI. Obtained first moments and spin contents at Q 2 = 1 GeV 2 in 


the NLO analysis 


with fixed aq 












distribution 




a = 0.5 




a = 1.0 


a = 1.6 


Aq 




-0.077 




-0.057 


-0.051 


A 5 




0.550 




0.533 


0.294 


AS 




0.123 




0.241 


0.276 
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TABLE VII. Obtained first moments. 



d i stxi "bn t/i on 


O 2 


LO 


NLO-1 


NLO-2 




1 HpV 2 
-L vje v 












n Q9fi ffivprH 
u.yzu ^nxeuj 


D Q9fi (ftvpr\~) 


fl Q9fi ffivprH 


Ad 




— D ^41 ('fivprl'l 


— D ^41 ('fivpH'l 
u.o^ii 1 11-x.eu. i 


— D ^41 ('fivpH'l 


An 




— Dfi4 

U.UU^i 


— n dsq 

u.uoy 


— n 0^7 


A ^, 




U.ool 


U.DoZ 


U.ooo 


o p 
in 




0.144 


110 


128 


n n 
9l 




— D 067 

U.UU i 


—0 069 

U. UUi/ 


—0 051 

U. U(J ± 


9i 




n 036 


01 Q 

U. U ± 


036 




O VJC V 








An 








n Q3n 

u. you 


Arl 




— *344 


— 343 

U.04:0 


—0 344 


An 




U.UU / 


U.UOj 




^9 




1 Ql/j 

l.ol4 


U.ouo 


U.yzU 


9i 




U. 14:0 


u. ± ±o 


n 1 37 


9i 




— n 068 

u. uuo 


—0 075 

U.U I (J 


—0 056 


a d 
9i 




n 035 

u. uou 


090 


038 




10 HpV 2 








An 




D Q94 

U. yZi^i 


U. i/OZi 


n Q3i 

U. cfO± 


Arl 




— ^4^ 


— 343 


— 34^ 


An 
L\q 




u.uuo 


u.uoy 


— D^IQ 
\j.\JOcf 


A 5 




1.524 


0.999 


1.077 


9i 




0.143 


0.120 


0.139 


9 n i 




-0.068 


-0.076 


-0.057 


gf 




0.035 


0.021 


0.038 
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FIGURES 
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FIG. 1. Q 2 evolution results for the proton structure function gf. The initial LO and NLO-1 
gi structure functions are evolved to those at 60 GeV 2 by the LO and NLO DGLAP evolution 
equations. 
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FIG. 2. Calculated LO and NLO spin asymmetries A\ for the proton are compared with the 
experimental results by the SMC, SLAC-E143, and HERMES collaborations at x « 0.117. The 
theoretical curves are obtained by using our LO and NLO-1 fitting results at Q 2 =l GeV 2 . 
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FIG. 3. Comparison of our calculations with the experimental asymmetry A\(x, Q 2 ) data for 
the (a) proton, (b) neutron, and (c) deuteron. Our results are obtained at Q 2 =5 GeV 2 with the 
optimum parameters in Tables II (LO) and III (NLO-1). The NLO and LO results are shown 
by the solid and dotted lines, respectively. 
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0.06 




(c) 

FIG. 4. Experimental spin-dependent structure functions xgi(x, Q 2 ) are compared with our 
LO results for the (a) proton, (b) neutron, and (c) deuteron. Our fitting results are calculated at 
Q 2 =l, 5, 20 GeV 2 by using the LO evolution equations with the optimum parameters in Table 
II, and they are shown by the dashed, solid, and dotted curves, respectively. The experimental 
data are obtained from the A\(x,Q 2 ) data and the F2(x,Q 2 ) calculated with the unpolarized 
GRV distributions and R\ggo{x, Q 2 ). 
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0.03t 

- - 1 GeV 2 * E142 
0.02 — 5 GeV 2 a £154 NLO 



20 GeV 2 D HERMES 

0.01 




-0.03^ ■ ■ — ■ ■ — 

0.001 0.01 0.1 

X 



(b) 
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FIG. 5. Experimental data of xg\(x,Q 2 ) are compared with our NLO-1 results for the (a) 
proton, (b) neutron, and (c) deuteron. The notations are the same as those in Fig. 4. 
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X 
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• SMC 



0.1 



■a 
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0.5 



-0.5 



-1.5 



0.001 



- Ad v (x)/d y (x) 
SMC 



0.01 




0.1 



(b) 

FIG. 7. Our LO ratios Au v /u v and Ad v /d v are compared with the SMC data in (a) and (b), 
respectively. Our results are calculated at Q 2 = 10 GeV 2 by using the polarized distributions in 
the LO. 
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0.001 



FIG. 8. The antiquark distributions of transformed SMC ("SMC") and LSS(1999) are com- 
pared with our NLO-1 distribution at Q 2 =l GeV 2 . 
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X . 
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FIG. 9. The x m i n dependence of AS(i m j n ) = f 1 AYl(x)dx is compared with the recent 
parametrizations of SMC and LSS(1999) at Q 2 =l GeV 2 . 
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